Plant Communications

The combination of target capture sequencing (TCS) with low-coverage whole genome sequencing (lcWGS), an approach known as Hyb-Seq, has allowed the integration of natural history collections into the genomics revolution, transforming biodiversity research. To implement Hyb-Seq, a collection of genomic targets, often nuclear orthologs, is needed to design probes for TCS. In flowering plants, the universal Angiosperms353 probe set has been proven resolutive at multiple evolutionary scales, with caveats. Malpighiales is known to be one of the most challenging flowering plant orders to resolve. Within this order, the clusioid clade ([~]2.2K species, 94 genera, five families) is no exception. To resolve phylogenetic relationships in this recalcitrant clade, we design a custom probe set, the Clusioids626 kit, composed of 39,936 120-mer probes targeting 626 nuclear orthologs ([~]6.6M nucleotides). This probe set includes all Angiosperms353 targets and 273 clusioid-specific ones, carefully chosen taking copy-number, length evenness, and phylo-informativeness into account. We test our probe set on 70 accessions representing all families and tribes in the clusioid clade. On average, 50.4% of TCS reads mapped to our targets, recovering a median of [~]600 orthologs. Relationships for all clusioid families are fully resolved for our nuclear targets. Additionally, 105 plastid coding DNA sequences were retrieved from the lcWGS fraction. A strong cyto-nuclear conflict was detected. The Clusioids626 kit performs better than the universal Angiosperms353 enrichment panel alone. Our kit design workflow can be extended into other lineages for which a universal probe set exists but more resolution is needed.

Poplar seed fibers cause environmental and health concerns, yet their developmental mechanisms remain poorly understood. Here, we constructed a high-resolution spatiotemporal transcriptomic atlas of female poplar capsules by integrating single-nucleus and spatial transcriptomics. We delineated the developmental trajectory of seed fibers, confirming their origin from placenta cells, and identified three functionally distinct fiber cell subtypes involved in initiation, metabolic support, and elongation. Weighted gene co-expression network analysis (WGCNA) identified several hub transcription factors, including PtoMYB, PtoHDT1, PtoEIF6 and PtoPDF2, that may serve as key regulators of fiber development. Our study provides a cellular-resolution framework for understanding trichome development in woody perennials and offers candidate targets for functional characterization toward breeding low-fluff poplar cultivars. HighlightsO_LIA spatiotemporal transcriptomic atlas of poplar capsule development is constructed at single-cell resolution C_LIO_LIFiber cells originate from placenta cells and comprise three functionally distinct subtypes C_LIO_LIProvides molecular targets for breeding low-fluff poplar cultivars to mitigate environmental pollution C_LI

Sacha inchi (Plukenetia volubilis L.) is an emerging woody oilseed crop prized for its high alpha-linolenic acid (ALA) content. Despite its nutritional and economic value, the lack of high-quality genomic resources has hindered genetic improvement and the elucidation of its unique polyunsaturated fatty acid and lipid biosynthetic pathways. In this study, we report a high-quality, chromosome-scale genome assembly of sacha inchi with a total length of 710.62 Mb, integrated from Illumina, PacBio, and chromosome conformation capture (Hi-C) technology. The genome harbors 37,570 protein-coding genes, and 379.86 Mb (53.45%) of repetitive sequences. Phylogenomic analysis reveals that sacha inchi diverged from its closest relative Ricinus communis, [~] approximately 36.2 million years ago. Comparative genomics indicates that sacha inchi experienced only ancient whole genome duplication events. To elucidate the mechanisms governing ALA biosynthesis and triacylglycerol (TAG) accumulation in sacha inchi seeds, we performed temporal transcriptome profiling across six seed development stages. Our findings demonstrate that high TAG content is primarily driven by the sustained expression of biosynthetic genes and low activity of degradation genes during mid-to-late seed development. Notably, while genes encoding stearoyl-ACP desaturases (SADs) maintain the precursor pool, the expression of genes encoding fatty-acid desaturase 2 (FAD2) and fatty-acid desaturase 3 (FAD3) is positively correlated with the final accumulation of C18:2 and C18:3 fatty acids. We also identified lncRNAs as potential epigenetic regulators of these key pathways. This high-quality genome provides a critical foundation for elucidating the molecular mechanisms of seed growth and development in sacha inchi.

Fitness costs of plant disease defence are often subtle and difficult to quantify. In this study, we therefore used comparative high-throughput phenotyping in two independent facilities to assess growth, morphology and physiology of potato (cv. Desiree) with high time-resolution monitoring different defence mechanisms under pathogen-free conditions. Plants were either treated weekly with the resistance inducers {beta}-aminobutyric acid (BABA; 10 mM) or potassium phosphite (KPhi; 36 mM) or comprised six transgenic lines expressing late blight resistance genes (single Rpi genes or a three-gene stack) or reduced jasmonate perception (StCOI1-RNAi). Over four weeks, image-derived traits revealed consistent cross-facility effects for plant height and colour: BABA treatment increased plant height but reduced canopy area and induced a paler greenness signature, whereas KPhi caused minimal and transient growth effects. Chlorophyll fluorescence at the NaPPI facility indicated reduced vitality (Rfd_Lss) in BABA-treated plants and increased Rfd_Lss following KPhi, while maximum PSII efficiency was largely unchanged. Several transgenic lines showed somewhat reduced above-ground biomass. Enzyme activity profiling produced distinct treatment and genotype signatures, but was strongly modulated by facility conditions that overrode these specificities. Overall, high-throughput phenotyping robustly detected subtle growth-defence trade-offs across platforms. HighlightHigh-throughput optical phenotyping validated across two independent research facilities reveals that stacked resistance genes and resistance inducers in potato trigger subtle growth trade-offs. Graphical abstracts O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=97 SRC="FIGDIR/small/713143v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@89df47org.highwire.dtl.DTLVardef@1a1ce64org.highwire.dtl.DTLVardef@1f52f0dorg.highwire.dtl.DTLVardef@1e41c35_HPS_FORMAT_FIGEXP M_FIG C_FIG Experimental timeline for high-throughput plant phenotyping platforms. Created in BioRender. Poque, S. (2026) https://BioRender.com/nmkve7g

Bioavailability of iron, an essential micronutrient to plants, is low in alkaline or calcareous soils, which are prevalent across semi-arid production regions. Breeding efforts to increase tolerance to iron deficiency chlorosis (IDC) in sorghum, a major crop of semi-arid regions, are confounded by spatial variation of stress severity in field trials. Here we developed and validated two high-throughput phenotyping approaches to address this challenge, with multi-spectral aerial imaging in the field and a controlled-environment assay to isolate the effects of iron bioavailability. In the field, severity and uniformity of stress are highly predictive of genetic signals for IDC tolerance (R2 > 0.6 for soil pH metrics and H2). Plot-level data filtering for stress conditions based on control genotypes successfully addresses field spatial variation (unfiltered H2 = 0.18 vs. filtered H2 = 0.4). The controlled-environment assay proxies field stress using iron sources with differential bioavailability, evidenced by high heritability ( H2 = 0.98) and phenotypic differential for hybrid control genotypes that matches field performance. Finally, we show that assay phenotypes are suitable for genome-wide association studies in global germplasm. Together, these field and lab phenomic approaches can be deployed to understand genetics of IDC tolerance and develop crops resilient to alkaline soils. HIGHLIGHTStress severity and uniformity greatly impact detection of genetic signals underlying iron deficiency chlorosis tolerance in sorghum. A controlled-environment assay reduces spatial heterogeneity and improves assessment of tolerance genetics.

Hi-C data is commonly used for reference-free de novo scaffolding. However, with the rapid increase in high-quality reference genomes, reference-guided workflows are now more practical for assembling large numbers of target genomes without relying on costly and labor-intensive Hi-C sequencing. Recently, a pangenome graph-based haplotype sampling algorithm was introduced to generate personalized graphs for target genomes. Such graphs have strong potential as references for reference-guided contig scaffolding. Here, we present noHiC, a reference-guided scaffolding pipeline supporting key steps of plant contig scaffolding. A distinctive feature of noHiC is the nohic-refpick script, generating a best-fit synthetic reference (synref) from a pangenome graph that is genetically close to the target contigs. This enables the integration of genetic information from many references (up to 48 in our tests) without using them separately during scaffolding. Synrefs showed advantages over highly contiguous conventional references in reducing false contig breaking during reference-based correction. Additionally, nohic-refpick can be combined with fast scaffolders (ntJoin) to rapidly produce highly contiguous assemblies using synrefs derived from pangenome graphs. The noHiC pipeline, used alone or in combination with ntJoin, can generally produce assemblies that are structurally consistent with public Hi-C-based or manually curated genomes. The pipeline is publicly available at https://github.com/andyngh/noHiC. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=82 SRC="FIGDIR/small/712436v1_ufig1.gif" ALT="Figure 1"> View larger version (9K): org.highwire.dtl.DTLVardef@40bd8forg.highwire.dtl.DTLVardef@5d2bbborg.highwire.dtl.DTLVardef@e214a3org.highwire.dtl.DTLVardef@b90b06_HPS_FORMAT_FIGEXP M_FIG C_FIG

Common buckwheat (Fagopyrum esculentum Moench) is a globally cultivated pseudocereal with a high nutritional quality and economic value. Due to its self-incompatibility, common buckwheat exhibits a high level of heterozygosity, making genome assembly challenging. Consequently, reference-level haplotype-resolved assemblies of common buckwheat are scarce, hindering research and genomics-assisted breeding. Here, we present a near-complete, chromosome-level, haplotype-resolved assembly of a common buckwheat F1 genotype (named Tuka), generated using a trio-binning approach that integrated parental Illumina short-read data with PacBio HiFi and Hi-C data from Tuka. The Tuka assembly comprises two haplomes, Tuka_h1 and Tuka_h2, both showing high contiguity (contig N50 of 76.68 Mb and 84.57 Mb, respectively), high completeness (assembly sizes of 1.28 Gb and 1.23 Gb with BUSCO scores of 96.9% and 96.8%, respectively), high base-level accuracy (QV of 59.08 and 63.03, respectively), and few gaps (35 and 30, respectively). This near-complete assembly of Tuka serves as a valuable genomic resource for common buckwheat, enabling advanced genomic analyses and accelerating research and breeding using state-of-the-art genomic tools.

Sesame (Sesamum indicum) is one of the earliest domesticated oilseed crops and is valued for antioxidant lignans that stabilize oil quality. However, the genomic and evolutionary history of the genus Sesamum, including the origin of its allotetraploid relative S. radiatum and the diversification of lignan metabolism, remains poorly understood owing to limited chromosome-scale genomic resources. Here we present chromosome-level genome assemblies for three wild Sesamum species, two Ceratotheca species and a Japanese sesame cultivar to reconstruct genome and karyotype evolution across the Sesamum-Ceratotheca complex. Comparative analyses show that the derived x=16 lineage originated from an ancestral x=13 karyotype through chromosome fission, fusion and translocation, whereas another x=13 lineage underwent extensive restructuring associated with retrotransposon expansion. Phylogenomics places Ceratotheca within the x=16 Sesamum clade and reveals that S. radiatum originated through hybridization involving a C. sesamoides-like ancestor. The antioxidative lignan gene CYP92B14 was reintroduced via the BB progenitor, linking hybridization with restoration of oil-stabilizing metabolism during sesame evolution.

Hundreds of proteins in the cell require iron (Fe) or Fe-containing cofactors to function. However, how Fe2+ or Fe3+ are specifically allocated to each of these proteins in plant cells remains largely unknown. It has been proposed that Fe metalation could be driven by specific interactions with Fe-shuttling proteins known as Fe-chaperones. Here, we present the first family of plant Fe2+-chaperones (ICHAPs) with orthologues in dicots and monocots. The role of these proteins in Fe distribution to Fe-dependent metabolic processes has been illustrated using symbiotic nitrogen fixation in Medicago truncatula root nodules. ICHAP1 is a soluble Fe2+-binding protein that interacts with plasma membrane Fe2+ transporter NRAMP1, but not with symbiosome Fe2+-transporters. ICHAP1 mutants present altered Fe distribution in cells and they cannot fix nitrogen. A second family member, ICHAP2 is required to target Fe2+ to symbiosomes, as it accepts Fe2+ from ICHAP1 and interacts with symbiosome Fe2+-importer VTL8, but not with NRAMP1. These results indicate a path for Fe2+ allocation from the plasma membrane to the symbiosome through specific protein-protein interactions and Fe2+ exchange from NRAMP1 to ICHAP1, to ICHAP2, and to VTL8.

Iron (Fe) deficiency is a widespread disorder limiting global soybean (Glycine max (L.) Merr.) production. Although root exudation is a key adaptive mechanism for Fe scarcity in species like Arabidopsis, a detailed chemical characterization of soybean exudates is lacking. Here, we examined the accumulation and secretion of phenolic compounds in soybean roots and their correlation with intraspecific tolerance to Fe-deficiency chlorosis. Seven soybean genotypes with contrasting tolerance, derived from U.S. breeding programs, were analyzed. Root exudates from Fe-deficient soybean plants solubilized ferric oxide. We identified and quantified 28 coumarin-type phenolics, with catechol methylsideretin as the predominant component. Although the qualitative coumarin profile was consistent across all genotypes, Fe-efficient lines secreted these compounds at higher levels or earlier during Fe deficiency than Fe-inefficient lines. The efficient genotype A7 showed coordinated upregulation of coumarin biosynthesis and secretion, whereas this response was weaker in the Fe-inefficient genotype IsoClark. Catechol methylsideretin concentrations strongly correlated with the ability of root exudates to mobilize Fe from ferric oxide. The conserved phenolic profile, together with divergence from those reported in non-legume species, suggests lineage-specific adaptations and ecological roles beyond Fe mobilization. These results highlight genotype-dependent exudation as a determinant of soybean Fe-deficiency tolerance, with implications for breeding. HIGHLIGHTIron deficiency induces soybean root exudates containing predominantly catechol methylsideretin which mobilize iron; genotypes differing in Fe efficiency show conserved qualitative but contrasting quantitative coumarin profiles.

AO_SCPLOWBSTRACTC_SCPLOWNitrogen fertilizers are essential for crop productivity but cause environmental harm, necessitating the development of cultivars that thrive under limited nitrogen. This study investigates the transcriptomic response to nitrate in Arabidopsis thaliana (a model dicot), Brachypodium distachyon (a model Pooideae), and Hordeum vulgare (barley, a domesticated Pooideae) to identify conserved and species-specific molecular mechanisms. Using RNA-seq after 1.5 and 3 hours of nitrate treatment, we found that core nitrate-responsive biological processes - such as nitrate transport, assimilation, carbon metabolism, and hormone signaling - are largely conserved across species. However, comparative analysis at gene level based on orthology revealed specificities between the species. For instance, rRNA processing was uniquely stimulated in Arabidopsis, while cysteine biosynthesis from serine and gibberellin biosynthesis were specifically regulated in Brachypodium and barley. Orthologs of key nitrate-responsive genes (e.g., NRT, NLP, TCP20) exhibited variable regulation, reflecting potential adaptations linked to domestication or nutrient acquisition strategies. These findings highlight the importance of integrating model and crop species to uncover targets for improving nitrogen use efficiency in cereals. The study provides a pipeline integrating gene ontology and orthology analyses to compare transcriptomic responses between species.

Climate change increasingly threatens global Capsicum (pepper) production. Accelerating the deployment of climate-resilient cultivars requires effective use of genetic diversity conserved in genebanks. We implement a "turbocharging" strategy in Capsicum by integrating genome-wide association studies and genomic prediction in a core collection (n = 423), followed by genomic prediction across the global collection (n = 10,250) using the core as a training population. We generated genomic estimated breeding values (GEBVs) for 31 high-accuracy traits (r > 0.5) encompassing hyperspectral phenotypes (heat/control), agronomic performance (heat/control) and fruit quality. To enhance accessibility and decision-making, we developed a large language model (LLM) integrated application that enables flexible, preference-based selection of candidates. By narrowing the parental decision space, this framework streamlines screening of large germplasm collections while balancing climate resilience, quality attributes and market demands. Our approach provides a scalable decision-support system to accelerate climate-resilient Capsicum breeding and maximize global genetic resources.

Pepper (Capsicum annuum L.) is a recalcitrant species regarding shoot regeneration, a trait that serves as a major bottleneck for the application of genetic engineering tools. In this study, comparative genetic analysis between a rare high-regeneration cultivar and a common low-efficiency cultivar identified a single nucleotide polymorphism (SNP) in PHYTOCHROME A SIGNAL TRANSDUCTION 1 (CaPAT1) that determines shoot regeneration efficiency. The T478C SNP in the high-efficiency cultivar converts a stop codon into an Arg codon, leading to translational read-through into the neighboring gene and forming an intact GRAS domain. This SNP-mediated formation of full-length CaPAT1 is essential for its dimerization. Notably, the overexpression of CaPAT1T478C in multiple low-efficiency cultivars, including both hot and bell peppers, significantly improved both shoot regeneration and transformation efficiency in the transformed T0 generation. These findings demonstrate the pivotal role of CaPAT1 in enhancing shoot regeneration and provide a robust strategy to overcome recalcitrance in pepper.

Plant amino acids function as both pathogen nutrients and essential drivers of systemic immunity. The regulation of amino acid homeostasis through transporters is a essential for mounting a robust and coordinated immune response in plants during pathogen infection. Using systems biology and integrative network science, we investigated bacterial virulence in Arabidopsis. By comparing gene coexpression networks of effector-triggered susceptibility (ETS) and pattern-triggered immunity (PTI), we uncovered a plant amino acid-related processes specifically linked to ETS. Integrating time-series transcriptomics, protein-DNA interactions, and mathematical simulations, we identified ANAC046 as a transcriptional regulator of amino acid processes during ETS. Single-cell RNA-Seq revealed that amino acid transporters are primarily expressed in companion and mesophyll cells, while functional validation confirmed ANAC046s roles in promoting susceptibility. Further integration of transcriptome and interactome data showed that amino acid-related genes interact with key immune hub proteins. Network topology analysis enabled the characterization of seven additional genes involved in plant defense. To support community-wide research, we developed MIData, an open-access platform for pre-analyzed Arabidopsis networks. Together, our findings demonstrate the power of systems-level approaches in uncovering hierarchical regulatory mechanisms underlying plant susceptibility to bacterial pathogens.

The narrow genetic base of cultivated tomato (Solanum lycopersicum L.) represents a major constraint on crop improvement, necessitating the exploitation of wild relatives to broaden allelic diversity. Here we present SABER (Solanum lycopersicum Allele Biodiversity Enriched Resources), a novel eight-founder Multiparent Advanced Generation Intercross (MAGIC) population that, for the first time, incorporates the Galapagos wild relative Solanum cheesmaniae as a founder alongside seven elite S. lycopersicum lines. Following a structured crossing scheme and Single Seed Descent advancement, F6 recombinant inbred lines were genotyped at 5,850 high-confidence SNP markers using Single Primer Enrichment Technology (SPET). Population structure analyses confirmed low residual heterozygosity, limited substructure among offspring, and successful introgression of S. cheesmaniae alleles across all twelve chromosomes. Mapping performance was validated through three Mendelian traits with known genetic determinants, all of which resolved to genomic positions consistent with the literature. QTL mapping for quantitative agronomic traits identified known loci for fruit epicarp and flesh color, and two novel QTL for days to flowering, number of leaves before flowering, and soluble solids content. Together, these results demonstrate that SABER is a powerful and reliable platform for high-resolution QTL mapping and candidate gene discovery, and establish a replicable framework for integrating wild germplasm into multiparental tomato breeding resources

Concerns over climate change have intensified the demand for stress resistant crops like hybrid wheatgrass (HWG; Elymus hoffmannii, StStStStHH), a perennial forage species known for its exceptional salt and drought tolerance. However, hexaploidy and high heterozygosity have complicated efforts to resolve its genomic structure and evolutionary history. Here, we present high-quality, haplotype-resolved, chromosome-level genome assemblies for HWG (CDC Saltking) and its putative progenitor, bluebunch wheatgrass (Pseudoroegneria spicata). By integrating PacBio HiFi and ultra-long Oxford Nanopore sequencing with Hi-C scaffolding, we assembled the 10.7 Gb HWG genome into 21 pseudochromosomes per haplotype. Our phylogenomic analysis redefines the origin of the H subgenome, positioning it as an intermediate between Old-World Hordeum marinum (sea barley) and Hordeum brevisubulatum. Notably, we identified significant chromosomal rearrangements, including a unique duplication on St chromosome 4. Transcriptome analysis across multiple tissues revealed a pronounced expression dominance of the H subgenome. This dominance was not associated with reduced LTR density, suggesting that selective pressures for rapid adaptation of the latest subgenome entrant may drive its dominance. Finally, using the f-branch statistic, population genomic analysis of 189 accessions representing eight Elymus and Pseudoroegneria species revealed extensive reticulate evolutionary relationships and identified P. spicata as a major, asymmetric genetic donor within the wheatgrass complex. These resources provide a foundational framework for future genomic research and genetic improvement in grasses and for the introgression of stress-tolerance traits into cereal crops such as wheat. Key MessagesDevelopment of world-first high-quality chromosomal-level haplotype-resolved genome assemblies of hexaploid HWG and diploid progenitor, Pseudoroegneria spicata, enabled the identification of the subgenome origins. This study resolved the evolutionary placement of the St genome and clarified the history of polyploidization and hybridization in HWG. Homeolog expression bias in the H subgenome likely reflects selective pressure favoring greater gene retention and upregulation of functionally important genes, thereby enhancing hybrid fitness. Population structure analysis distinctly differentiates P. spicata, E. repens, E. hoffmannii from other European Pseudoroegneria species. The findings reveal the complex patterns of interspecific gene flow and population dynamics within the Elymus and Pseudoroegneria species.

Drought is a major constraint on the productivity of durum wheat across Mediterranean and North African regions. To elucidate the mechanisms underlying drought resilience, we employed a combination of scenario-controlled phenomics and flag leaf transcriptomics across ten durum wheat genotypes. These included the Tunisian landraces Chili and Mahmoudi, seven breeding lines, and the reference cultivar Svevo. The plants were grown to maturity under well-watered or long-term drought conditions in pots and rhizotrons, enabling a comprehensive assessment of growth, yield components, root architecture, physiological traits, and reaction norm plasticity. Drought markedly reduced performance, yet Chili and Mahmoudi consistently maintained superior biomass, grain number and intrinsic water use efficiency (iWUE). This was supported by balanced C/N allocation, strong osmotic adjustment, and the ability to sustain robust root systems under stress, albeit through partly divergent physiological strategies. Transcriptomic profiling revealed highly genotype specific responses, with drought tolerance unrelated to the number of differentially expressed genes. Instead, the landraces displayed distinct regulatory programs involving mainly photosynthesis protection, ABA-related transporters, osmotic adjustment pathways, and stress-responsive transcription factors. These mechanistic insights identify actionable physiological and molecular determinants of drought plasticity and provide high value targets for accelerating the breeding of climate resilient durum wheat. HighlightsIntegrated phenomics and transcriptomics revealed landrace-specific physiological and molecular mechanisms enabling superior drought resilience and identifying actionable targets for durum wheat improvement.

Strigolactones (SLs) are plant hormones regulating shoot branching and symbiotic interactions, but their trace-level abundance limits research and applications. Here, we optimized a Nicotiana benthamiana transient expression system for SL production by tuning agroinfiltration parameters and co-expressing rate-limiting carotenoid biosynthetic genes. Overexpression of Zea mays PSY1 or an Arabidopsis PSY-GGPS11 fusion increased carlactone production over 2-fold and enhanced downstream SL accumulation. Using this platform, we discovered that sorghum cytochrome P450 SbCYP728B35 catalyzes conversion of 5-deoxystrigol to sorgolactone, revealing a previously unknown function. These results establish metabolic engineering of precursor supply as an effective strategy for boosting SL production and demonstrate N. benthamiana as a robust system for pathway elucidation and biotechnological synthesis of bioactive strigolactones.

Computational tools, including biological language models (LMs), show substantial promise in predicting the impact of genetic variants on plant fitness. However, validating variant effect predictions (VEP) requires experimental populations where genetic variation consists of discrete point mutations rather than segregating recombination blocks. In this study, we generated a novel population of Brachypodium distachyon mutant lines to evaluate the accuracy of VEP at single-base resolution. These lines were advanced through single-seed descent for five generations (M1 to M5), with whole-genome sequencing performed at M2 and M5 and phenotypic measurements recorded at M3 and M4. Using state-of-the-art VEP models, we predicted the functional impact of missense protein-coding variants and gene-proximal non-coding variants. We validated these predictions by estimating the effect of mutations on whole-plant measurements (burden tests) and their probability of fixation from M2 to M5 (purging tests). Among missense variants, the protein LM ESM showed superior predictive accuracy compared to the bioinformatic standard SIFT and the genomic LM PlantCAD. Notably, the relationship between VEP scores and allele fixation suggested a log-linear relationship between VEP scores and variant fitness. Among gene-proximal variants, PlantCAD appeared more accurate than supervised models of regulatory activity, such as chromatin accessibility (a2z) and RNA abundance (PhytoExpr). Collectively, our findings highlight the utility of state-of-the-art VEP tools as predictors of fitness and demonstrate the potential of mutant populations to evaluate computational tools for precision breeding applications.

Salt stress alters plant development, including the floral transition, but regulation of timing of flowering by salt is poorly understood at the molecular level. To identify genetic loci regulating the floral transition under high soil salinity, we performed a genome-wide association study (GWAS) in Arabidopsis thaliana and identified natural variation at the UGT74E1-UGT74E2-BT3 (UUB) locus that correlates with bolting time specifically in response to salt stress. Genetic analysis revealed BT3 as a novel repressor of the floral transition in control conditions. Similarly, the putative IBA glycosylases UGT74E1 & UGT74E2 delay the floral transition in control conditions. Furthermore, we identified that IBA homeostasis regulators TOB1 and ECH2/IBR10 play a key role in the floral transition, and that ECH2/IBR10 are required for the early flowering phenotype of the ugt74e1/ugt74e2 double mutant, indicating that UGT74E1 & UGT74E2 delay flowering by altering IBA homeostasis. A pangenome analysis of the UUB locus revealed variation in the occurrence of the DNA transposon SAUERKRAUT (SKRT). CRISPR-mediated SKRT deletion in Col-0 affected gene expression both within and outside the UUB locus and caused a salt-dependent delayed floral transition. The delayed bolting phenotype of the skrt-2 mutant also depends on ECH2/IBR10 function, indicating that SKRT accelerates the floral transition by altering IBA homeostasis. Finally, targeted demethylation of SKRT resulted in delayed floral transition under salt stress. Taken together, our data show a role for SKRT and its DNA methylation levels in the salt-dependent bolting time response in Arabidopsis, revealing a novel molecular mechanism to control flowering in adverse conditions.